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Abstract 

The asymptotic behaviour near phase transitions can be suitably characterized by the scaling of 
As/Q 2 with e = 1 — T/T c , where As is the excess entropy and Q is the order parameter. As As is 
obtained by integration of the experimental excess specific heat of the transition Ac, it displays little 
experimental noise so that the curve log(As/Q 2 ) versus loge is better constrained than, say, log Ac 
versus loge. The behaviour of As/Q 2 for different universality classes is presented and compared. 
In all cases, it clearly deviates from being a constant. The determination of this function can then 
be an effective method to distinguish asymptotic critical behaviour. For comparison, experimental 
data for three very different systems, Rb2CoF4, Rb2ZnCl4 and SrTiOs, are analysed under this 
approach. In SrTi03, the function As/Q 2 does not deviate within experimental resolution from 
a straight line so that, although Q can be fitted with a non mean-field exponent, the data can 
be explained by a classical Landau mean-field behaviour. In contrast, the behaviour of As/Q 2 
for the antiferromagnetic transition in Rb2CoF4 and the normal- incommensurate phase transition 
in Rb2ZnCl4 is fully consistent with the asymptotic critical behaviour of the universality class 
corresponding to each case. This analysis supports, therefore, the claim that incommensurate phase 
transitions in general, and the A2BX4 compounds in particular, in contrast with most structural 
phase transitions, have critical regions large enough to be observable. 
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1. INTRODUCTION 



The observation of fluctuation-driven critical phenomena in structural phase transitions 
is subject of long and persistent controversies. Although power law behaviours near phase 
transitions can be "universal" , the temperature range where this asymptotic behaviour can 
be observed within experimental resolution (the so-called asymptotic or critical region) is 
not. It depends in general on the particular features of each system and can even vary 
for different physical quantities. It is no wonder, therefore, that the size of the critical 
region in structural phase transitions has been a matter of discussion for several decades. In 
general, one should expect it to be smaller than in the pure order-disorder systems, the spin 
hamiltonians, that represent the corresponding universality classes. However, the Ginzburg 
criterion has been equally invoked to predict typical critical regions of size too small to 
be of any experimental significance, 0] as to justify critical regions of the order of 0.1 in 
reduced temperature (T/T c ) [2f (an interval that would greatly surpass that of the actual 
spin-hamiltonian representatives of some of the universality classes) . In addition, structural 
defects may produce strong deviations from the universal asymptotic laws predicted for 
ideally perfect crystals. 0, 0,0] 

Authors expecting insignificantly small critical regions have argued repeatedly that the 
experimental observation of power laws around the phase transition with non-mean-field 
exponents may have a rather simple explanation, e.g. in terms of the Landau theory when 
sixth-order or higher order terms are taken into account in the Landau expansion |lj. For 
instance, a 2-4-6 Landau potential: 

AG = 1-AeQ 2 + ]bQ 4 + \c(f (1) 
2 4 o 

with temperature independent parameters A,B,C and e = 1 — T/T c yields indeed a mean 
field critical behaviour for the order parameter (Q oc e 1 / 2 ) with a mean-field exponent 
(3 = 1/2, but the asymptotic behaviour is limited to e < B 2 /AAC. This temperature 
interval can be beyond any observation if the fourth order coefficient is small enough, a 
typical situation happening if the system is close to a tricritical point. In e-intervals of the 
order of 0.01-0.1, or even larger, the behaviour of the order parameter can then be fitted to 
power laws with "effective" exponents which are intermediate between the mean field value 
1/2 and the one corresponding to mean field tricritical behaviour (1/4). 

More generally, it became clear that any experimental observation of expected asymp- 
totic critical behaviour for the order parameter requires further confirmation through cross- 
examination of its consistency with the behaviour of other physical parameters near the 
transition point. The correlation of the order parameter with the transition excess entropy 
has been recently used for this purpose [6]. The excess entropy As is expected to be propor- 
tional to the square of the order parameter Q not only in the mean field asymptotic regime 
(P = 1/2, a = 0), but also for any behaviour described by a Landau potential like equa- 
tion The relation would be true as long as the terms higher than quadratic in the 
Landau potential are temperature independent and the quadratic term has a linear tem- 
perature dependence. This particular relation between entropy and order parameter is in 
principle broken in the asymptotic regime for all other universality classes. Its experimen- 
tal check can be then an effective way to distinguish between genuine fluctuation-driven 
critical behaviour and pure Landau behaviour. Indeed, this proportionality relation has 
recently been shown for SrTiC>3, in contradiction with a common belief that this compound 
exhibits in a significant temperature interval asymptotic critical behaviour corresponding to 
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the Heisenberg-3D or the Ising-3D universality class |6|, LZL ISl] - The analysis of As, instead of 
checking directly the expected critical behaviour of the excess specific heat Ac, has the fun- 
damental advantage of using an integrated quantity with a smooth temperature dependence 
that is independent of experimental statistical errors in the measurement of Ac. 

In the present paper, the approach of comparing As and Q 2 for analysing the thermal 
behaviour at phase transitions is further investigated and applied to several materials. The 
expected behaviour of As/Q 2 for the different universality classes and the Landau case, 
is discussed and compared with that derived from experimental data for Rb2CoF 4 (anti- 
ferromagnetic phase transition), R^ZnCLi (normal- incommensurate phase transition) and 
SrTiOa (antiferrodistortive transition). 



2. ENTROPY-ORDER PARAMETER RELATION IN THE CRITICAL REGION 

The asymptotic entropy-order parameter relation for the different universality classes 
follows from the critical laws for specific heat and order parameter. It can be expressed as 



Q 

with 

h 

and 



2 oc e K As oc Q 2K (2) 




2(3 a > f . 

a < [6) 
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In order to get these expressions the normal asymptotic behaviour of the order parameter, 
characterized by the exponent j3, is assumed while that of the excess specific heat is supposed 
to be: 

Ac oc (5) 

a 

Note that the specific heat is not reduced to a pure power law as it includes the cases of 
negative a values and the logarithmic divergence for a — > 0. [ol flfj| 

The resulting exponents k and k' for the different universality classes and the Landau 
behaviour are listed in table |U Non-classical exponents « and k' differ significantly from 
the mean-field/Landau trivial values associated with the above mentioned proportionality 
law. Nevertheless, these deviations may be rather difficult to detect in a typical As vs. Q 2 
plot. In figure H we present these plots for the different universality classes assuming that, 
in an e interval of 0.1 and within experimental resolution, both the order parameter and the 
excess specific heat follow their asymptotic critical behaviour (see equation (j5J) for that of 
the excess specific heat). 

The excess entropy in figure Q has been calculated by numerical integration of the ex- 
pression: 

A^/^dT-/^' (6) 
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TABLE I: The critical exponents for several universal classes. The dimensionality of the "inter- 
action" is labeled d. The dimensionality of order parameter is labeled n. Exponent a is that of 
the specific heat, 8 is the order parameter exponent, In the limit e — > 0, k and k' are given by 
equation (j3J) and equation K eT stands for the "effective" exponent obtained in an e-interval of 
10 -3 — 10~ 2 for the theoretical assymptotic behaviour of the models. 
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FIG. 1: Theoretical behaviour of the excess entropy vs the square of the order parameter for various 
universality classes. Landau mean field theory gives a straight line passing through origin. Order 
parameter is normalized to unity at e = 1, while entropies have been renormalized so as to match 
at Q 2 = 0.2. Upper labels stand for e in the XY-3D universality class, they are also approximately 
valid for Heisenberg-3D and Ising-3D but not for Ising-2D or for Landau behaviour. 

this integration leads generally to smooth curves As(T) even if Ac contains significant 
statistical noise. 

Classical Landau behaviour is represented in figure Q] by a straight line through the 
origin. The Ising-3D, Heisenberg-3D and XY-3D classes yield curves different but close to 
a straight line. Therefore, in practice, the deviation from the Landau behaviour is difficult 
to assess quantitatively through these plots. Only in the case of the Ising-2D class the curve, 
being a fourth order parabola, clearly deviates from the proportionality law. The different 
relationships between As and Q 2 become more obvious in a log-log plot of As/Q 2 vs e. The 
plots corresponding to the curves of figure^are presented in figureEl The deviation from the 
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0.001 0.01 0.1 
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FIG. 2: As/Q 2 versus e in a log-log plot for the same data used for curves of figure [I] The 
magnitude represented in y-axis is normalized to unity at e = 0.1 except that of Ising-2D which 
has been arbitrarily shifted upwards for the sake of clarity. 

horizontal Landau behaviour is then obvious for all universality classes. The larger deviation 
is seen for the Ising-2D case, while all other universality classes lead again to more similar 
curves. It also becomes evident in the figure that although the curves have been obtained 
assuming perfect asymptotic behaviour for Ac and Q up to e = 0.1, the magnitude As/Q 2 
is far from being in the asymptotic regime represented by the power law of equation (jHJ) even 
at e values so low as 0.0005. The curves deviate rather smoothly from straight lines and in 
fact their approximate slopes are far from the asymptotic values given by equation (jSJ) and 
listed in tabled Indeed, the "critical" region for the excess entropy is much smaller than for 
Ac, since the integration represented by equation © yields the sum of several competing 
terms with similar but different powers of e. Hence, according to figure El the asymptotic 
behaviour for As/Q 2 represented by equation (|2j) is not experimentally accessible in any 
case. Nevertheless, the curves in figure El can be considered "universal" in the sense that 
they should be satisfied for any system so long as its order parameter and specific heat are 
following the asymptotic critical behaviour of the corresponding universality class. In the e 
interval [0.001-0.01] the curves can be approximated to straight lines, their slopes correspond 
then to "universal" effective exponents n e for the approximate law As/Q 2 oc e Ke r restricted 
to this temperature interval. The exponents k £t obtained from the linear fit of the curves 
in figure El in the mentioned temperature interval are listed in table |U 

3. EXPERIMENTAL CASES 

We consider now under this approach three very different systems: strontium ti- 
tanate (SrTiOs), rubidium tetrafruorocobaltate (Rb2CoF 4 ) and rubidium tetrachlorozincate 
(Rb 2 ZnCl 4 ). 
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FIG. 3: Excess entropy versus excess intensity of X-ray difraction experiments around the SrTi03 
phase transition at about 106 K. In the upper axis approximate values for reduced temperature 
are indicated 

3.1. SrTiO, 



The antiferrodistortive phase transition of SrTiOa has been widely considered one of 
the few examples of a structural phase transition with a clear observable critical region. 
Miiller et a/.Q could fit their EPR data assuming a power law for the order parameter 
with (3 = 0.35 in an e interval of the order of 0.1. This behaviour was then identified as 
the expected asymptotic behaviour of the Heisenberg-3D class. Recently, the excess entropy 
was determined and compared with the order parameter data; a mean field linear correlation 
was observed j<§], thus the critical region might be beyond observation. For comparison with 
earlier specific heat data see Ref. LLJ. 

For the sake of comparison with the two following systems, the As-Q 2 relation obtained 
in Ref. for SrTiC>3 is reproduced in figure El in a format similar to next figures. 



3.2. Rb 2 CoF 4 



This compound exhibits an antiferromagnetic phase transition at 107 K. This system is 
believed to belong to the Ising-2d universality class. The structure of Rb2CoF 4 is isomorphic 
to that of K2MF4 which is a well known two-dimensional anti-ferromagnetic system. The 
two-dimensional character is related to the presence of coupled anti-ferromagnetic planes 



separated by planes with weak interplanar effective interactions. Specific heat [12j and 
the reduced sublattice magnetization (order parameter) obtained from neutron scattering 
experiments k| have been reported to be consistent with the critical behaviour predicted for 
this universality class. In figure 0] we present the correlation of the excess entropy entropy 
As and the square of the order parameter (the reduced sublattice magnetization a) for this 
material. The entropy has been calculated upon integration of the specific heat measured 
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FIG. 4: The excess entropy of the antiferromagnetic phase transition in RD2C0F4 versus the square 
of the reduced magnetization, as obtained using data from Samuelsen(l^| and Ikeda et a/.0] Upper 
labels represent reduced temperature. 

by Ikeda et a/.Jl2|. The magnetization was determined from neutron scattering experiments 
by Samuelsen|l3l]. The figure shows clearly how the entropy fails to be proportional to the 
square of the order parameter. 



3.3. Rb 2 ZnCl 4 

Rh^ZnCU presents an orthorhombic (Pnma) phased) above T c ~ 305 K. At this tem- 
perature it transforms through a continuous phase transition into an incommensurately 
modulated crystal with a modulation wave vector close to one third of a*. Further below, 
at 195 K the modulation wave vector locks into the exact value 1/3 a*, and the incom- 
mensurate phase transforms into an orthorrhombic commensurate ferroelectric phase with 
a three-fold unit cell along the modulation. This compound belongs to the family of type 
A 2 BX 4 , where this incommensurate instability is rather common, and has been extensively 
studied in order to test the theoretical predictions on incommensurate phase transitions. In 
1978 Bruce and Cowley[15j] argued that this transition should belong to the XY-3D uni- 
versality class. Since then, many workers have related their results to what is expected in 
the critical zone for this universality class. For instance, Walisch et al.^^ accurately de- 
termined from quadrupole perturbed NMR experiments the critical exponent of the order 
parameter with full agreement with the theoretical value for the XY-3D class. According 
to these results the critical zone would reach up to e values of the order of 0.07. Diffraction 
synchrotron measurements (l7| confirmed these results and a critical region of the order of 
0.1. Specific heat was also reported [18f to agree expected critical behaviour in a range up 
to e ~ 0.2. More recently, Haga et al.[l£$ made a more comprehensive study of the specific 
heat and fitted their data to an expression which includes corrections up to second order of 
the asymptotic law. The fitted e interval reached in this case up to e = 0.1. The existence 
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of observable large critical regions in incommensurate phase transitions, in comparison with 
normal structural phase transitions, was predicted and justified. |2(], |2l[ However, other au- 
thors maintain that there is no theoretical ground for expecting in incommensurate systems 
larger critical regions than in other structural phase transitions 0, H3| and have interpreted 
for instance other NMR experiments in this compounds within a Landau approach 23]. 

For the present study, a new independent measurement of the specific heat of F^ZnCU 
was performed in a conduction calorimeter. The calorimeter has been previously described 
*24l | and has been used successfully for the measurements of thermal properties of single- 



crystals 25[. A single crystal of 1.078 g (2.850 mmol) was measured. The sample was grown 
from aqueous solution and was annealed at the high symmetry phase and then cooled down 
at a low rate of 0.4Kh _1 while the heat capacity of the sample was recorded. The heat 



capacity agrees previously reported data [Hi |2y, |22j and they show a A-shaped peak at 
about 304 K. 

The background contribution to the specific heat was determined by fitting a smooth 
polynomial function to data above T c and data well below the lock-in phase transition located 
at Tj = 195 K. By subtracting the background contribution, the excess specific heat was 
determined and the excess entropy versus temperature was obtained through integration. 
This latter could then be compared with the temperature evolution of the width Av of the 
NMR frequency distribution measured by Walisch et al. [3| , which is a quantity proportional 
to the order parameter. Due to the sharpness of the specific heat anomaly we have assumed 
that T c corresponds to the temperature of maximum specific heat. The temperature origin 
of the NMR data was slightly modified to make T c coincide. Excess entropy was set to 
zero at e = 0. Finally, an interpolation function was used to get the values of the NMR 
frequency distribution width Au at the temperatures values were entropy was determined. 
Figure 03 presents the correlation between both quantities in a As versus Au 2 graph. The 
deviation from a mean field linear behaviour is rather small, but as shown in figure ^ XY 
3D asymptotic behaviour implies in this representation small changes with respect to the 
Landau mean field behaviour. 



4. DISCUSSION 

4.1. Comparative analysis of experimental data 

Figure IH1 depicts the log-log plot of As/Q 2 versus e for the three compounds discussed in 
sectional One can clearly see the deviation of both RbgZnC^ and Rb2CoF4 from the Landau 
mean field behaviour, in contrast with the case of SrTiC>3. In accordance with the simulations 
in section the curves in general can be approximated to straight lines when reduced to 
the interval of smaller e values. Their slopes yield "experimental effective exponents" K e , 
listed in table |Hj that can be compared with the theoretical ones listed in table |l] It can be 
stated that in the case of Rb2ZnCL the experimental effective exponent n e value coincides 
with the theoretical value of the XY-3D class considering its standard deviation. It should 
be noted that the values of similar for XY-3D and Ising-3D universality classes. 

However, the value (3 that fits the order parameter behaviour is consistent with only the 
XY-3D class. The K e value for Rb2CoF4 is close to k £t expected for Ising-2D class, with a 
deviation of the order of 5%; while in SrTiC>3 not only the experimental n e clearly disagrees 
with the theoretical value of the Heisenberg-3D class, but its zero value indicates a classical 
Landau behaviour. It should be stressed that a mere deviation from the linearity between 
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FIG. 5: Excess entropy versus square of the NMR frequency distribution width Av proportional to 
the order parameter in Rh^ZnCU. The straight line represents the Landau mean field behaviour. 
Labels at upper x-axis represents e. 

Crystal e fitting range Points n e 

SrTiO-3 lO^-KT 1 22 0.006(5) 

Rb 2 CoF 4 10~ 3 -10- 2 12 0.633(6) 

Rb 2 ZnCl 4 10~ 3 -10- 2 16 0.169(6) 

TABLE II: Experimental effective exponent K e , deduced from the log-log fit in figure El for the 
systems under study. Values of experimental n e should be compared to those expected theoretically 
for the universality classes K eT listed in table [I] 

As and Q 2 should not be taken as a proof of fluctuation-driven critical behaviour, since a 
generalized Landau potential with temperature dependent quartic and higher order terms 
could explain such behaviour. It is the quantitative agreement of this deviation (measured 
by the effective exponent n e ) with that expected for the corresponding universality class 
(k 6t ), together with the agreement of the fitted exponent (3 that proofs the consistency of 
an interpretation in terms of asymptotic critical behaviour for both the order parameter and 
the specific heat. 

4.2. Crossover in Rb2ZnCi4 

The curve of Rb2ZnCl4 in figure |H1 strongly deviates from its low e linear behaviour for e > 
0.01 and changes then rapidly into a near horizontal straight line. This suggests a crossover 
around this temperature from the XY-3D critical behaviour to a mean-field/Landau regime. 
However, this interpretation would be somehow in contradiction with the fact that the 
observed order parameter behaviour, used in fact for the calculation of this curve, does not 
suffer any kind of crossover at e ~ 0.01 and follows the j3 power law with excellent agreement 
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FIG. 6: Log-log plot of the ratio As/Q 2 — in arbitrary units — as a function of e for Rl^ZnCU, 
SrTiOa, and RD2C0F4. Straight lines show the linear fitting done in the intervals indicated in 
Table 2. Resulting effective exponents summarized in table ILT1 

up to much larger e values, of the order of 0.1. In any case, the curve clearly deviates for 
e > 0.1 from the theoretical curve for the XY 3D class in figure Eland this can then only 
be caused by the deviation beyond this e value of the experimental Ac from the theoretical 
asymptotic behaviour given by equation (JHJ) and used to produce figure El The reason for 
such deviation could be intrinsic, i.e. the critical region for the specific heat would be then 
of the order 0.01, while for the order parameter is about 0.1 in reduced temperature, but 
it is also possible that the critical region were sample dependent. In order to get a more 
comprehensive analysis of this phase transition we have also analyzed the measurement done 
by Haga et a/. [HI where the excess specific heat was fitted over 30 K, i.e. 0.1 in reduced 
temperature, using a complex expression involving corrections up to second order of the 
pure asymptotic regime: 



where a = -0.007, A = 8.153, D 1 = 0.040, D 2 = -0.026, B c = 8.201, A 1 = 0.524; where 
A, Di, D2, B c are expresed in arbitrary units. This curve is compared in figure with that 
obtained from a mere critical law as given in equation (J5|) . which essentially corresponds to 
the neglect of the coefficients D% and D2 in equation (J7|). 

The difference between both curves become obvious beyond e = 0.01, showing that for 
lower temperatures corrections to the asymptotic behaviour are significant. For comparison, 
we have computed again the curve of As/Q 2 using as specific heat data the values given by 
the function equation (JJJ) as proposed by Haga. The results are presented in figure |H1 

The agreement between both calculations is excellent. The approximate linear behaviour 
for e values smaller than 0.01 is reproduced and its slope gives a value of around 0.15 similar 
to that obtained from our specific heat analysis. On the other hand, a small shift along the 
e-axis of the crossover region can be observed and the behaviour for larger e values slightly 
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FIG. 7: The solid line represents the curve (6) used by Haga et a/.(l9(to fit the excess specific heat 
of Rt^ZnCU. The dashed line is a simple critical law as given by equation © with a = —0.007 in 
accordance with the XY-3D class. 

differs. These differences are most probably due to the choice of the baseline for obtaining 
the "excess" specific heat. The influence of the baseline is slight at small values of e but 
becomes very important when e increases. The results presented for RL^ZnCL in figure El 
correspond to a baseline described by a polynomial of minimal order fitting the specific data 
well above and below the phase transition. This polynomial can be significantly changed, 
by shifting the low-temperature fitting region to higher o lower temperatures, while keeping 
a physically consistent baseline. The effect of this change of baseline is shown in figure |H1 
It is slight for low e values, where the approximate linear behaviour remains unchanged, 
while the crossover to what one might consider a Landau regime (the horizontal part of the 
curve in figure IBJ) becomes less pronounced. As shown in the figure, the change is such that 
the calculated curve coincides with the theoretical one of the XY-3D class in a much larger 
e interval. Hence, a large uncertainty exists with respect to the actual e-interval where a 
crossover to the Landau behaviour does take place. But, the behaviour of the As/Q 2 (e) 
relationship for e < 0.01 confirming a XY-3D critical regime is robust with respect to 
modifications of the baseline. 
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